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Abstract. The Fermi Large Area Telescope (LAT) has revolutionized high-energy (HE) 
astronomy, and is making enormous contributions particularly to 7-ray pulsar science. 
As a result of the many new pulsar discoveries, the 7-ray pulsar population is now 
approaching 100. Some very famous millisecond pulsars (MSPs) have also been detected: 
J1939-h2134 (B1937-h21), the first MSP ever discovered, as well as J1959-h2048 (B1957-K20), 
the first black widow pulsar system. These, along with other MSPs such as PSR J0034— 0534 
and J2214-f3000, are rare among the pulsar population in that they exhibit nearly phase- 
aligned radio and 7-ray light curves (LCs). Traditionally, pulsar LCs have been modelled using 
standard HE models in conjunction with low-altitude conal beam radio models. However, a 
different approach is needed to account for phase-aligned LCs. We explored two scenarios: one 
where both the radio and 7-ray emission originate in the outer magnetosphere, and one where 
the emission comes from near the polar caps (PCs) on the stellar surface. We find best-fit LCs 
using a Markov chain Monte Carlo (MCMC) technique for the first class of models. The first 
scenario seems to be somewhat preferred, as is also hinted at by the radio polarization data. 
This implies that the phase-aligned LCs are possibly of caustic origin produced in the outer 
magnetosphere, in contrast to the usual lower-altitude conal beam radio models. Lastly, we 
constrain the emission altitudes with typical uncertainties of ~ 10% of the light cylinder radius. 
The modelled pulsars are members of a third 7-ray MSP subclass, in addition to two others 
with non-aligned radio and 7-ray LCs. 



1. Introduction 

For quite some time, there were only a handful of pulsars detected in the 7-ray waveband [T] . This 
small sample already provided a starting point for spectral and light curve (LC) modelling as 
well as preliminary population studies. However, the launch of the Fermi Large Area Telescope 
(LAT) [2] heralded an exciting new era, particularly for high-energy (HE) pulsar physics. 



1.1. J -ray pulsars galore 

Following the first Fermi Catalog which included about 1 500 7-ray sources [3j, the second 
Fermi Catalog is in now production, containing nearly 1 900 7-ray sources [U [5] including 
almost 100 7-ray pulsars [Si [7]. About a third of these are millisecond pulsars (MSPs), some 
having been found by observing non-variable, unassociated Fermi sources at high latitudes with 



radio telescopes and searching for pulsed radio signals [SI O [10] . Some of these MSPs exhibit 
the unusual phenomenon that the peaks of their 7-ray and radio LCs occur at the same observer 
phase ('longitude' 4>), i.e., they are phase-aligned. 

1.2. Traditional models 

Two classes of models have been invoked to explain HE pulsar radiation. The first class, low- 
altitude polar cap (PC) models [IIl[l2], assume that primary electrons are accelerated above 
the neutron star surface, and that magnetic pair production of curvature radiation or inverse- 
Compton-scattered 7-rays occurs in the intense B-fields close to the stellar surface. In the 
case that pair creation is suppressed along the last open field lines, which mark the boundary 
between the corotating closed field line zone and the active open one, a slot gap (SG) [13] 
may form, corresponding to a two-pole caustic (TPC) geometry [H] which extends from the 
stellar surface up to near the light cylinder. The second class of models are the outer gap (OG) 
models [15\ I16j . which assume the production of HE radiation along the last open field lines 
above the so-called null charge surface, where the Goldreich-Julian charge density becomes zero. 
In both classes of outer-magnetospheric models (SG and OG), the HE pulse profiles are the 
result of the formation of caustics - the accumulation of photons in narrow phase bands due to 
a combination of special relativistic effects and magnetic field line curvature. Also, the narrow 
gaps in these models require abundant production of electron-positron pairs which will be able 
to screen the accelerating E-field. 

For the case of MSPs, LCs and spectra have been modelled [l7l [18l [19] using a pair-starved 
polar cap (PSPC) model ^20j. The latter is similar to the traditional PC model, but in this case 
the number of pairs is not sufficient to screen the accelerating E-field, such that the accelerating 
region above the star is pair-starved. MSP spectra and energetics have also been modelled 
using an OG model |2H 122] . while LC modelling using an annular gap model can furthermore 
reproduce the salient features of the 7-ray LCs [23J. 

The outer magnetospheric models usually invoke relatively low-altitude 'core' and 'conal' 
radio beams centred on the magnetic axis, in addition to the HE radiation being produced in 
extended regions reaching up to the light cylinder. The difference in the location of the 7-ray 
and radio emission regions thus implies that the corresponding pulse profiles will have non-zero 
lags (phase offsets) between the 7-ray and radio LCs. 

1.3. Phase-aligned radio and ^-ray light curves: a new MSP subclass 

The LCs of the first 8 Fermi-detected 7-ray MSPs [H] have been modelled [25] , yielding two 
distinct MSP subclasses: those whose LCs are well fit by a standard OG or TPC model, and 
those whose LCs are well fit by a PSPC model. These fits are mutually exclusive. Importantly, 
this unexpectedly implied that there is copious pair production taking place even in MSP 
magnetospheres with their characteristic low B-fields. 

A third class of MSPs emerged with the discovery of PSR J0034— 0534 [26], which exhibited 
(nearly) phase-aligned radio and 7-ray LCs. Prior to this detection, such behaviour has only been 
observed for the Crab pulsar [27]. This phenomenon has now also been seen for PSR J1939+2134, 
PSR J1959+2048 [28], and PSR J2214+3000 [TO]. Some other MSPs, including PSR B1821-24 
and PSR J0737— 3039A, may also have phase-aligned radio and 7-ray LCs, on the basis that 
their radio and X-ray profiles are phase-aligned [23 132] • 

2. Geometric models predicting phase-aligned light curves 

The non-zero radio-to-7 lags of MSPs with phase-aligned LCs disqualify the application of 
standard models, since the phase alignment implies co- located emission regions. We consider 
two model scenarios, one where the 7-rays and radio both occur in extended, high-altitude 
regions, and one where they originate at low altitudes. (For more details, see [31j.) 



2.1. High- altitude models 

We use the same framework as previously [25], except that the radio emission region is now 
extended in altitude, and not assumed to be coming from a radio cone at a single altitiude. We 
furthermore limit the minimum and maximum radii of the radio and 7-ray emission regions, 
and assume constant emissivity as a function of altitude. These models are therefore called 
altitude-limited OG (alOG) and TPC (alTPC) models. Technically, the traditional OG model 
is a specific instance of an alOG model, and similar for the traditional TPC and alTPC models. 

2.2. Low- altitude models 

The so-called low-altitude Slot Gap (laSG) models represent an alternative possibility assuming 
a non-caustic origin of the emission. These are actually very-low-altitude geometric SG models 
resembling a hollow-cone beam close to the stellar surface. We modulate the emissivity as 
motivated by detailed radiation models, where the emissivity rises and falls exponentially along 
the B-field lines, peaking at a distance of about one or two stellar radii. Additionally, we also 
investigated the case where the emissivity is constant. 

3. Finding optimal LC solutions in multidimensional phase space 

The models described above involve multiple free parameters describing the emission region 
extent and width, as well as the pulsar geometry (inclination and observer angles a and Q). In 
order to pick statistically the best-fit parameters when comparing to the Fermi LG data, we 
have developed a Markov chain Monte Carlo (MCMC) maximum likelihood procedure |32j and 
applied this to the alOG and alTPC models for three MSPs as indicated below (see Table [TJ. 
The 7-ray LCs are fit using Poisson likelihood while the radio LCs are fit using a statistic. 
These two values are then combined. Using this method, we can also derive uncertainties on the 
minimum and maximum altitudes of the emission regions. 

4. Results 

4-.1. The alOG and alTPC models 

Increasing the lower limit of the emission region's altitude while keeping the upper limit fixed 
has the effect of growing the PC size (i.e., deleting a ring of radiation around the PC in the 
phaseplot, which indicates relative intensity per solid angle vs. C, and 4>). This suppresses 
low-level off-peak emission, especially for the alTPC case, so that the peaks become sharper. 
Conversely, decreasing the maximum radius while keeping the lower one fixed constrains the 
emission to a ring-like structure around the PC. As a result, the relative peak heights change in 
the corresponding LCs, and low-level features may become more prominent. 

Figure [2 shows an example of best alOG and alTPC LC fits for PSR J1939-h2134. The 
best-fit values for a and for the three MSPs we have modelled, are summarized in Table [U 



Table 1. Best-fit pulsar geometries for different models. 
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Figure 1. LC fits for PSR J1939+2134 using alTPC and alOG models. Panel (a) shows the 
7-ray data [28], while panel (b) shows the radio data. 

having typical errors of ~ 5 — 10°. Best-fit emission altitudes and gap widths are presented 
elsewhere [31]. 

Jf..2. The laSG models 

It is interesting that in the laSG models the leading peak is wider while the trailing peak is 
sharper, even for the low-altitude emission considered here. This means that the caustic effects 
start to appear, given the large corotation speed of MSPs already near their surface. By including 
emission from higher altitudes, the peaks become broader. Also, the peak phase separation and 
width may be altered by choosing different values for the minimum and maximum emission 
altitudes. The constant-emissivity assumption leads to block-shaped LCs, and is therefore not 
considered viable. For the modulated-emissivity case, the peak widths may be fine-tuned using 
fading parameters, while the peak separation depends sensitively on C. This class of models 
shows less variation in profile shape than the high-altitude ones, and there are also multiple 
combinations of the free parameters that give very similar LCs, so a 'best solution' is probably 
not unique. The different best-fit solutions usually are a trade-off between best fits for 7-ray vs. 
radio LCs. For this reason, we indicate two similar laSG solutions (labeled 'laSGl' and 'laSG2', 
differing only in fading parameters which set the emissivity fading properties) in Table (H 

5. Discussion 

Pulsars with radio and 7-ray LCs which are non-aligned in phase may be modelled by 7-ray 
emission regions extended over a large range of altitudes (OG / TPC / PSPC models) in 
conjunction with a conal radio beam at relatively lower emission altitudes. However, phase- 
aligned LCs require co-located 7-ray and radio emission regions. We studied the viability of 
reproducing such phase-aligned LCs using both high-altitude and low-altitude geometric models. 
We found that both classes of models could produce LC fits that can capture the most prominent 
features of the pulse profiles. 



5.1. Model constraints 

We have found best fits for the two angles describing the pulsar geometry, a and as indicated 
in Table [TJ These have typical errors of ^ 5 — 10°. In addition, we were able to limit the emission 
geometry for the high- altitude models, and found that the typical errors on the best-fit emission 
altitudes are ~ 10% of the light cylinder radius. Furthermore, the maximum radio emission 
altitude seems to be better constrained than the 7-ray one. Since the low-altitude models can 
produce very similar LCs for similar values of its free parameters, we regard the inferred values 
of these free parameters when applying the model to the data as less robust than the case of 
the altitude-limited models. 

5.2. Which model is preferred? 

While we used the MCMC technique to search for the optimal set of parameters that produced 
the closest match to the LC data for the case of the altitude-limited models, we have not yet 
implemented this technique for the laSG models. Since the parameter space has not been fully 
explored for the latter class of models, we cannot really claim to have found the best fit for 
the alSG models. It is therefore difficult to quantitatively favour one class of models above 
the other. However, we did calculate the likelihood of the manually-selected best-fit laSG LCs 
as well as the MCMC-selected best-fit LCs of the altitude-limited models, and found that the 
altitude-limited models give better fits than the laSG ones, while the alTPC models provide 
slightly better LC fits than the alOG models. This implies that the phase-aligned 7-ray and 
radio LCs are most probably of caustic origin, produced in the outer magnetosphere, and the 
radio emission is most likely originating near the light cylinder (see also [33j). 

5.3. Does causticity imply phase-alignment? 

By studying the visibility of two samples of 7-ray and radio pulsars, it was concluded [34J that 
radio and 7-ray beams must have comparable sky coverage, especially for pulsars with large spin- 
down luminosities. This implies that the radio emission should originate in very wide beams 
at a significant fraction of the light cylinder, motivating studies of high-altitude caustic radio 
emission. However, this scenario is limited in application to pulsars having nearly phase-aligned 
LCs (with phase differences of up to ^ 0.2), unless the radio emission region is significantly offset 
from the 7-ray one. The bulk of the 7-ray pulsar population exhibits quite large (~ 0.1 — 0.5) 
radio-to- 7-ray phase lags, so that caustic radio emission is probably is not ubiquitous in young 
pulsars. Furthermore, the 7-ray profiles are usually double-peaked, while the radio ones are 
single-peaked. Unless the radio emissivity has a strong altitudinal or azimuthal dependence, 
both radio and 7-ray LCs would be mostly double-peaked if they result from caustic emission, 
which is not observed. On the other hand, radio caustics may be more common for short-period 
MSPs, as there are many more examples of MSPs with phase-aligned LCs. Detailed population 
studies involving both young and old pulsars will provide more quantitative answers to this 
question. 

5.4. The potential of polarization 

The rotating vector model [35] generally provides a good description of radio polarization data 
(position angle as function of phase). However, this model assumes a static dipole B-field with no 
rotation effects included, and is probably not valid for fast-spinning MSPs which are expected to 
have significant B-field distortion. Caustic emission models predict rapid position angle swings 
with phase, coupled with depolarization [36], since the emission from a large range of altitudes 
and B-field orientations is compressed into a narrow phase interval to form the peaks. Rapid 
changes in the polarization angle and low levels of linear and circular polarization near the peak 
phases have indeed been seen for the MSPs modelled in this paper, and may be indicative of 
caustic effects. Polarization signatures will therefore be important to help discriminate between 



LCs produced by caustic emission (such as occurs in alOG / alTPC models) or not (e.g., in the 
laSG model). 

6. Conclusions 

Future studies should include development of full radiation models which will be able to 
reproduce both the multiwavelength LC shapes, polarization properties, as well as the energy- 
dependent behaviour of the spectra of the 7-ray pulsars. Ways in which to induce abundant 
pair creation in low-B MSP magnetospheres, such as using offset-dipole j^lj or multipole [21] 
B-fields, will need to be investigated. 
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